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Independent studies of attenuation of X-band (8.5 GHz) and Ka-band (32 GHz) 
radio signals when traversing Martian dust were carried out by the authors. These analy- 
ses turned out remarkably similar. The computational method is essentially that of 
T. S. Chu but uses observed optical depth at 0.67 microns rather than “visibility” as 
the measure of optical attenuation from which to derive the microwave attenuation. 
An awkwardness in the approach is that the size distribution of Martian dust particles 
is not well known , but the mean is probably around 4 microns, whereas in the terres- 
trial case it is nearer 10 microns. As a consequence , there will be a larger tail of particles 
still in the Mie regime in the Martian case as compared to the terrestrial one. The com- 
putational error will , therefore, be somewhat larger for Martian than Earth-bound dust. 
Fortunately , the indicated attenuations are small enough for the worst case (1.3 dB at 
32 GHz) that the error is academic. 


I. Introduction 

This article makes use of the approaches and results of 
studies of terrestrial sand and dust storms for the rather 
analogous problem on Mars. Dust storms on Mars are, if any- 
thing, more impressive than on Earth. Typically, in spring in 
the southern hemisphere, dust storms will start, grow, move 
northward, and frequently encircle the globe. In 1971, there 
was a huge storm which obliterated the surface features of 
the entire planet. 
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translates to: L/L Q = 10 log 0.031 = 15.086 dB, which i 
usually shortened to 15 dB. An approach developed by Chu 
(Ref. 6), and used by Ansari and Evans (Ref. 1) and Gold- 
tiirsh (Ref. 8), incorporates the visibility V into the expres- 
sion for microwave attenuation coefficient A (2): 
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II. Terrestrial Approach 

In the case of the Earth, the most convenient yardstick for 
a dust storm is “visibility,” which is defined as the distance 
at which a mark disappears against the background. According 
to Middleton (Ref. 12), at the distance of disappearance, 


wher 


r = particle radius, m 
e = e r -je. = complex dielectric constant 
X = wavelength, m 
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As shown above, particles are assumed to be uniform in 
size and to have known permittivity. The development assumes 
single scatter. Jhis-is not bad #t microwave, but the optical 
'-visibility lifrft i^well into the multiple scatter regime. 

Particle size distributions for four dust storms in Khartoum 
have been published by Ghobrial (Ref. 7), who finds parti- 
cles with sizes ranging from 0.1 to 300 microns (0.0001 to 
0.3 mm). Wind tunnel measurements (Ref. 9) indicate that 
0.08 mm is the most likely size of particle to be picked up by 
wind from the Earth’s surface. 


III. Mars 

Mars has long captured the public fancy as possibly the 
most hospitable planet, next to Earth, for life. This is not so 
unreasonable, as can be seen in Table 1 . However, the tempera- 
ture is considerably colder, the atmospheric pressure at the 
surface less than 1% that of the Earth, and the atmosphere 
almost totally lacking in oxygen (95.3% C0 2 , 2.7% N 2 , 1.6% 
argon and perhaps 0.13% 0 2 ; Ref. 13). Mars’ orbit around the 
sun is considerably more eccentric than that of the Earth 
(206.7 X 10 6 km at perihelion and 249.1 X 10 6 km at aphe- 
lion) with the result that the insolation is some 45% greater 
during perihelion (late spring in the southern hemisphere) than 
at aphelion (late spring in the northern hemisphere). 

Mars has been visited by spacecraft from both the US and 
USSR as is seen in Table 2. The most recent visits were by the 
two US Viking missions, and the data presented here came 
largely from that source. 

Martian dust storms tend to originate in the late spring or 
summer in the southern hemisphere when the Martian surface 
temperatures and the temperature gradients are highest. The 
greatest storm on record started in September of 1971 and 
became planetwide (Ref. 11). Great storms (planetary scale) 
occur once or twice each Martian year (Ref. 19). Quite typi- 
cally, a storm will move northward and reach planet-encircling 
proportions and last several weeks. The dust has been reported 
as high as 50 km above the surface. Optical depths as great 
as 6 at 0.67 microns have been measured and a dust density 
scale height of 10 km proposed (Refs. 15 and 16). The size 
and permittivity of Martian dust remain controversial. Fig- 
ure 1 portrays a particle density distribution (Ref. 10) typical 
of those found in papers published in the late 1970s. The 
mean of the particle radius is 0.4 microns. Also shown is the 
geometric cross-section distribution, GN(r), which is more 
directly pertinent to extinction calculations. Recent workers 
have favored somewhat larger particles, around 2 microns 
(Ref. 19 and Footnote 1) or greater (Ref. 5). The favored 
substances for the dust are basalt and montmorillonitic clays. 


Tabulated dielectric constant at 50 microns (the longest wave- 
length available) for the latter is 2.18 -/0.14. 1 

Water was detected in the soils at the two Viking Lander 
sites, but the integrated water in the Martian atmosphere 
does not exceed 100 precipitable microns (Ref. 4). 


IV. Attenuation by Dust 

A simple slab model of a Martian dust storm is shown in 
Fig. 2. The assumption is made that the reflected components 
of the electric field E r and E tr are much less than the trans- 
mitted component E t (E t » E r E t » E tr ). We define a 
complex propagation constant^: 

K c = 0-/a (1) 

where 

2nn r 

j3 = & Q n r - -r — = phase constant (2) 

A o 

27 m. 

a = K Q n. = -r— ^ - = attenuation constant (3) 
A o 

and 

n. = n r -jn. = complex refractive index 



where 

X Q = free-space wavelength 

One can then write for the electric field (Ref. 2) 

E = if f (0, 1,0) exp {/(cor - Kz)} (4) 

where E f (0 y 1, 0) is the complex amplitude of the transmitted 
plane wave with electric vector polarized in the y direction and 
travelling downward (positive z direction). Substituting (1) 
into (4) yields 

E = E f (0,1 ,0) exp (-az) exp { j(cjt - 0z) } (5) 


*R. A. West, Jet Propulsion Laboratory, private communication, 
1986. 
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which may also be written 


E = E f (0,1 ,0) exp y rj exp j j(u>t - | (6) 

where r is the optical depth. The attenuation coefficient A in 
decibels is then given by: 

A (dB) = 8.686 a (nepers) (7) 

A (dB) = 8.686 K Q n. dB/m (8) 

Attenuation from dust = 4.343r dB (9) 

A. Chu’s Development 

We now turn to the work of Chu (Ref. 6) who wrote a 
much quoted article on “Effects of Sandstorms on Microwave 
Propagation.” Chu considered sand of radius 0.01 to 0.1 mm 
(10 to 100 microns) with dielectric constant 2.5 (1 - /0.01) 
to 10 (1 - /0.01). He invoked the Rayleigh approximation at 
microwave and the very large sphere approximation in the 
visual range. The attenuation coefficient a is defined as 


The question of multiple scatter needs to be raised. Van de 
Hulst (Ref. 18) suggests that whenever the optical depth 
exceeds 0.3 multiple scattering becomes a factor. In this 
case, the limit of visibility occurs at r = 3.47, well into the 
multiple scatter regime. However, this limit occurs at the 
point where a mark disappears into the background, and it 
can be argued that this is the coherent component that one is 
tracking and that single scatter still applies. 

The introduction of (12) is done by combining (11) and 
(12) and solving for N: 

a a 

N = Q“ fl3) 

6.5 [(4/3 )to 3 ] V ’ 

where a Q is the optical attenuation coefficient in dB/m. 

The next step is to make use of an expression presented, 
but not derived, by Van de Hulst (Ref. 18) for the effective 
refractive index m of a scattering medium: 

m = 1 - iS(0) 2nNK~ 3 (14) 


a 



N(a)C(a) da 


where 5(0) is the forward scattering function, and K is the free 

( 10 ) 

space phase constant. In the Rayleigh regime 5(0) is given by 


where 

C(a ) = extinction coefficient 
N(a) = number density 
a = particle radius 

For sand of radius a 


a = 4.34 


(4/3)ra 3 


x 3.2 5SQ . 

KO = y— dB / m 


(id 


where Q ext = C/ira 1 is the extinction efficiency, and S = (4/3) 
m 6 N is the fraction of the volume actually made up of sand. 
As the size of particles become large with respect to wave- 
length, Q ext becomes asymtotic to 2 (Refs. 17 and 18): 

Q ext = 2 ,*»A ( 12 ) 

As can be seen in Fig. 3, if the sphere is at all conducting, the 
extinction becomes very close to 2 by the time 2na/X = 15 or 
a = 2.4X. Hence, in the case treated by Chu, if we take 
0.6 microns for the wavelength of white light, any particle 
larger than 1.5 microns would meet the requirements of (12). 


s(0) ■ * 3 (ffl ) ‘*\ k> (ffj) 2 < 15 > 

The second term is negligible for X » a as is the case at 
radio frequencies. Substituting (13) and (15) into (14) and 
using the terminology of (1) through (9), the phase shift is 
given by 


= g V 




deg/m (16) 


and the attenuation coefficient becomes 


K(Im m) = a Q a 


Im 


(Kij 


(8.68) dB/m (17) 


B. Application to the Martian Problem 

For a uniform slab the total attenuation and total phase 
shift for the wave of relations (4) through (6) is obtained by 
multiplying (16) and (17) by z, the distance traversed. For 
total attenuation (17) becomes 


K(lm m)z = yj (% z)a 



(8.68) dB/m (18a) 
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and, substituting 4.34 r for a Q z, we may write 

K(lmm) = jjTa Im (ffy) (37.7) dB/m (18b) 

This relation (18b) now gives the microwave attenuation 
for uniform spheres of known radius a and dielectric constant 
e for a slab of constant particle density N and agrees with 
Eq. (13) of Goldhirsh (Ref. 8). Martian dust undoubtedly 
follows some particle distribution such as shown in Fig. 1, 
but the uncertainty as to the mean size is such that the as- 
sumption of a single particle radius is as good as any. Similarly, 
the dielectric constant can only be guessed at, but can rea- 
sonably be bounded. The dust density probably follows an 
experimental decay with height, so the equivalent slab will be 
used. 

C. Equivalent Slab for a Martian Dust Storm 

Pollack (Ref. 15) has suggested a scale height for dust 
particle density, H d , of 10 km. Assuming H d is a constant 
over the region of interest, we may write 

| ^0 exp (- -^jdh = N 0 H d (19) 

where the right-hand side is obtained by performing the indi- 
cated integration; and N 0 is the surface density of the dust. 
Thus, it is seen that the Martian dust storm may be repre- 
sented by a slab of height H d = 10 km and density N 0 . A 
spherical model of such a slab is shown in Fig. 4. 

D. Calculation of Microwave Attenuation 


For zenith angles of 0° to 80° we use: 

t = T z sec0, 6 <80° (21) 

For zenith angles greater than 80° we need to consider refrac- 
tion effects. The Martian atmosphere is 95% C0 2 ; C0 2 at 
0°C and 1 atmosphere and for white light has an index of 
refraction of 1.000450. In terms of refractivity N = {n - 1)1 0 6 ; 
this may be written 

^(C0 2 )= 121.3 p(mb)/JTK) (22) 

If p = 6 mb and r=280 K on the Martian surface, the 
surface refractivity is 2.6N units as compared to about 300 
on the surface of Earth, and refractive bending will be a 
second-order effect and can be ignored. Using the geometry 
in Fig. 4 for zenith angles 6 greater than 80° or elevation 
angles 0 between 0° and 10°, the optical depth r is 

T = /rjl^d + r o ) 2 - ('o cos 0) 2 ] 1/2 - ''o sin ^ |> 

for 0 = 90° -6 <10° (23) 

Shown in Fig. 5 is the zenith angle dependence for the 
“best guess” and “worst case” calculations using relations 
(20) and (23). Dust storms under the worst case scenario only 
become significant for ray paths nearly tangential to the 
surface of Mars. Shown in Fig. 5 is the dependence of atten- 
uation on choice of particle radius and dielectric constant. 

V. Discussion 


Performing the indicated complex arithmetic on (18b) 
leads to 


A = 54.62 


X (e, + 2) 2 + e] 


where e = e r -j€. and a, the particle radius, and X are in 
meters. 


Attenuation through Martian dust storms has been com- 
puted for 8.5 and 32 GHz one-way radio transmission for the 
maximum expected optical depth (r = 6). These are first- 
order calculations based on single scatter. Due to uncertain- 
ties in the Martian dust parameters - particle size distribution 
and dielectric constant - “best guess” and “worst case” 
situations have been selected. Attenuation, even at 32 GHz, 
will not be significant, except for rays almost tangential to 
the Martian surface. 
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Table 1. Some physical properties of Earth and Mars 



Property 

Earth 

Mars 

l. 

Mean distance from sun, km X 1G 6 

149.6 

227 9 ! 206 ' 7 
1 249.1 

2. 

Period of revolution, days 

365.28 

687 

3. 

Rotation period, hours 

23.9344 

24.6321 

4. 

Equatorial diameter, km 

12,756 

6,787 

5. 

Mass (Earth = 1) 

1 

0.108 

6. 

Density (Water =1) 

5.5 

3.9 

7. 

Atmosphere 

Nitrogen 

Carbon dioxide 



Oxygen 

Nitrogen, Argon 

8. 

Water vapor (integrated column), cm 2 

2.2 

0.001 

9. 

Water content in surface 

0-30% 

-1% 

10. 

Mean atmospheric pressure, mb 

1013.25 

6 

11. 

Mean surface temperature, °C 

22 

-23 

12. 

Mean particle size in dust storms (diam), 

0.5 to 100 

0.1 to 10 


microns 



13. 

Real part refractive index 

[2.5 + 0.5 (%H 2 0)] 


(clay 8.5 GHz) 



14. 

Imaginary part refractive index 

[0.02 + 

0.15 (%H 2 0)] 


(clay 8.5 GHz) 


15. 

Real part refractive index 

[2.5 + 

0.3 (%H 2 Q)1 


(clay 32 GHz) 


16. 

Imaginary part refractive index 

[0.06 + 0.1 (%H 2 0)] 


(clay 32 GHz) 



Notes: 1. Data for 1 to 8 from Ref. 14. 




2. Values for 13 to 16 derived for data for clay from Table 1, Ref. 1. 


Table 2. Spacecraft missions to Mars 8 


Mission 

Arrival date 

Type 

United States 



Mariner 4 

July 15, 1965 

Flyby 

Mariner 6 

July 30, 1969 

Flyby 

Mariner 7 

August 4, 1969 

Flyby 

Mariner 9 

November 14, 1971 

Orbiter 

Viking 1 

June 19, 1976 

Orbiter/ Lander 

Viking 2 

A n i Qiz 

/ j i. > ' v 

Orbiter ^ Lander 

Soviet 



Mars 2 

November 27, 1971 

Orbiter/ Lander 

Mars 3 

December 2, 1971 

Orbiter/ Lander 

Mars 4 

March 1974 

Orbiter 

Mars 5 

March 1974 

Orbiter 

Mars 6 

March 1^74 

Lander 

Mars 7 

March 1974 

Lander 

a From Ref. 3. 
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PARTICLE RADIUS, A 
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N(r) = ar a exp[- (a/r)(r/r m K ] 
for a = 2, y ~ 1 12, r m = 0.4 fj.cn 
(Ref. 17} 
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2 — A 
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= Q a + Q s 
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°a 

= EXTINCTION 
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n = 1 .29 

°s 

= EXTINCTION 
COEFFICIENT 
FOR SCATTERING 


PARTICLE RADIUS, 

Fig. 1 . Martian dust particle distribution 


Fig. 3. The extinction efficiency Q e for various types of spheres 
shown in the range 0 < x < 5 (Ref. 7) 


SLAB OF DENSITY Nn, HEIGHT H h = 10 km 


-d = 260.7 km 



tq = 3,393.5 km 



Fig. 4. Model of equivalent Martian atmosphere 


Fig. 2. Simple slab model of a Martian 
dust storm 









ONE-WAY ATTENUATION, dB 


* * 


V 



Fig. 5. Zenith angle dependence of dust attenuation in Martian 



r ARTICLE nAClUS, /Ain ( * 0 ^ mm) 


Fig. 6. Dependence of attenuation on choice of particle radius and 
dielectric constant 
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